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^ : 1 Introduction 
o 

^ ■ The two-fermion bound system is an attractive subject of atomic (positronium, hydrogen 
atom) and sub-atomic physics (deuteron, mesons). Despite these systems are rather simple 
the study of two-particle bound states is challenging and still remains a source of progress in 
^ , quantum theory. Last decade significant efforts were undertaken to describe such objects, 
^ I and progress has been achieved in both, solving the corresponding equations [l]-[5] and 
calculating the experimentally measured observables [6]-[T0]. 

The homogeneous Bethe-Salpeter (BS) equation nowadays remains a powerful tool to 
investigate the relativistic bound state problem. Recently the problem of rigorous treating 
the BS equations received renewed interest, and several successful attempts were made to 
^ I reconsider solving algorithms for the BS equation |Tl] and to make its reduction to Light 
Front form more transparent [4J. A good example of a consistent study of this subject can 

■ be found in Ref. [T2] . 

^ . In fact, our understanding of the mathematical properties of bound states within any 

^ ' relativistic approach is far from perfection. The BS equation itself is a quite complicated 
mathematical object, and the technical problem of solving it is still of principal value. In 
^ ■ the present paper we propose an efficient and promising method to solve the BS equation 

■ for fermions involving interaction kernels in the form of one boson exchange supplemented 
with the corresponding form factors. It is based on employing the basis of hyperspherical 
harmonics for expanding the partial amplitudes and kernels. We show that this new tech- 
nique allows one to utilize many advantages in understanding the BS approach. Basically, 
the current study is inspired by the results reported in [3]. We explore the structure of 0"*" 
and 1"^ bound states for different couplings studying in details the convergence of solutions 
and corresponding eigenvalues. In particular, on a basis of the introduced rigorous method 
to solve the BS equation it becomes possible to analyze in details the problem of stability 
of bound states within the BS approach. It is worth mentioning that this phenomenon 
was also considered in [3] in the framework of Light-Front Dynamics, and it is completely 
similar to the bound state collapse in non relativistic quantum mechanics for potentials 
behaving like ~ 1 /r^ . 



O 



2 Homogeneous Bethe-Salpeter equation for fermions 

In the present study, we generalize the method described in [HI [131 [H] to homogeneous BS 
equations for spinor particles. To provide a clear illustration, we present here the results of 
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solving the equation for the vertex function 

Gip) = ^ / y(p^ r(i) ^(^i) m s{h) f (2) (1) 

in the ^So channel with V{p, k) = g'^/[{p — k)"^ — /i^ + is] as one meson exchange interaction 
kernel for scalar, pseudoscalar and vector meson exchanges with corresponding vertices 
T{i),i = 1,2. The meaning of the introduced quantities p, k, ki, k2 is the following: fci 2 = 
P/2 ± k, k = {ko, k),p = {po,p) are the relative 4-momenta, and P = (M, 0) is the total 
4-momentum of the bound state in its c.m.s. Spinor propagators of constituent particles 
with equal masses m are 

^(fc)= ^^r, . , s{k)^cs{kfc 



^2 _ _)_ ^5 ' A;2 _ ^2 _|_ ' 

with C standing for charge conjugation matrix, C = ^7072, k = 'j^k^. In general, the 
vertices include the meson-nucleon form factor, which will be further considered as 

In ([1]) the BS vertex function Q{p) being the matrix 4x4 should be expanded over the 
proper set of matrices for the given channel. The standard choice for them is the set of 
the p-spin angular momentum vector harmonics Taip), where index a includes not only 

LSJ momenta but also p-spin quantum numbers, which are denoted by ++, , e and o 

|15j . For convenience in our further considerations we introduce another equivalent set of 
spin-angular matrices instead of p-spin basis. For ^5*0 channel the following functions may 
be chosen: 

ri{p) = \l,, TM = ^^ol„ Tab) = -^|j^7o75, ^^(p) = -^75- (3) 
This basis is orthonormal, i.e. 

dnpTr [r„(p)T„+(p)] = 



and the partial expansion can be written as 

G{po,p) = Y.9niPo, \p\) %{p). QniPo, \p\) = I dilpTr [g{p,,p)r+ {p)]. (4) 

After the partial expansion and upon performing the Wick rotation we can expand the 
partial vertex functions and the partial interaction kernels in hyperspherical functions 
using the formula [TT] 



1 1 ^ ~ 

VE{p,k) = - , ,.2 I ^2 = ^^Vy^^fe k)ZnlmiuJp)Z*i^{uJk), (5) 

\P ) E ' nlm 



where {k,p)E = k^p^ + {k,p), k = ykl + k"^ is the 4-dimensional absolute value, and 
"^fc = (x? ^5 0) are the angles of vector k = {k^, k) in 4-dimensional Euclidean space. The 
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Table 1: Pseudo-probabilities of partial components in the state with given M, i.e. their 
contributions to normalization condition. 

hyperspherical harmonics are Znimix^G^'P) = ^ni{x)'^im{0, 4>) (for details see e.g. [T^). 
Corresponding hyperspherical expansions for ^5*0 amplitudes are given by 

oo oo 

9i,2{^Pa,\p\) = J29i,2{p)X2j-2,o{Xp), 93{m,\p\) = J29iip)X2j-i,i{Xp), (6) 

oo 

94imAp\) = J29i{p)X2jAxp)- (7) 

By using these decompositions one can obtain the final system of 1-dimensional integral 
equations for the coefficient functions This system will be explicitly shown and dis- 

cussed in details separately [16J. What is important now, is that its numerical analysis does 
not require large computer resources. This set can be easily transformed to the system of 
linear equations. For this aim, firstly, the infinite summation over hyperspherical compo- 
nents should be limited to some finite value N^ax- Secondly, to calculate the integrals a 
reliable integration scheme is required. Applying Gaussian quadrature formula, one can get 
the system of linear equations with the sought functions defined in the mesh points [11] , 

X = XAX, (8) 

where X = g^, and the column 

= i[{9m}l'i]f=r, • • • , [{^7kfe)}£"i]f=r) 

represents the sought solution in the form of a group of sets of partial wave components 
gi^, a = 1, Nc^j = 1, N^ax specified on the integration mesh of order Nq- For 0"*" (1+) 
state we have Nc = 4: {Nc = 8). The matrix A is obtained as a product of partial kernels, 
the Jakobian of the transition to the new variables, the weights of the Gaussian mesh etc. 
The dimension of A is x A^, where A^ = NcNcNmax- 

3 Results 

In this short communication we are able to present only the most indicative results of 
the numerical treatment of the BS equations within our method. First of all, the set of 
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Figure 1: Dependence of the coupling con- 
stant a = g^/4:TT from the mass of the 
bound state M in different approaches. 



Figure 2: Functions gl,j = 1, 2, 3 for the 
cut-off A = 500 GeV/c and = 30 (sohd 
hne) and g^ = 48 (dashed hne). 



equations ([8]) has the solution only if det(AA — 1) = 0. This condition allows us to connect 
the coupling constant g^ and the mass of bound state M, i.e. for any given value of g^ the 
mass M can be calculated, and vice versa. The results of such calculation for a = g'^/^n are 
shown in Fig. [T]for the case of scalar meson exchange. The solid curve corresponds to our 
calculations in the BS approach, dashed curve represents the results obtained within Light 
Front Dynamics [3], dotted one - nonrelativistic calculations for the Yukawa potential. 

The spectrum of bound states obtained in this way demonstrates the customary non- 
relativistic features. Together with the ground states, the exited levels of the system can be 
found (for an alternative method of calculation see e.g. Ref. [IT]). Like in the nonrelativistic 
picture, solutions of the ground states (i.e. set of partial vertex functions gi, . . . ,g4) do not 
have nodes in |p| whereas the excited levels are described by vertex functions having zeroes. 

Besides, for any given mass of exchanged meson /x the bound state (at the ground level) 
can appear in the considered system only starting from some finite value of the coupling 
constant (7^j„. For example, at fi = 0.15 GeV we have g'^^J^ = 4.023, which corresponds to 
some minimal depth of the potential, where the bound state still exists. 

It is obvious, that for weakly bound states for the fixed binding energies B of order 
of a few MeV coupling constants are approximately equal. Thus, for the value B = 1 
MeV aBs = 0.362, and aiFD = 0.331. But for B of order of hundreds MeV coupling 
constants essentially differ. For example, for the value B = 100 MeV we have the ratio 
<^Bs/c(LFD ~ 1-3. The explanation of such a behaviour we found in the role of ^Sq~ 
component in the total BS solution. It is seen from the Table 1 that its contribution to 
the normalization condition is negligibly small for weakly bound states but increases very 
rapidly with the increasing B, in contrast to the contributions of other partial states. Thus, 
the role of the ^Sq component is repulsive, which leads to an increasing coupling constant 
for the same bound mass M in comparison with the nonrelativistic or LFD formalisms, 
where " " components are absent. 

Another important question to be touched upon is the problem of stability of the bound 
states within the BS approach. In the present context stability means the existence of the 
solutions in ([T]) without any cut off in the vertices. In its turn, such a solution exists only if 
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Figure 3: Masses of bound states in ^5*0 chan- Figure 4: Masses of bound states in ■^Si-^Di 
nel as functions of cut-off parameter A. channel as functions of cut-off parameter A. 

it does not depend on the parameters of calculation like Nmax, No, etc. In general, we found 
that convergence of our method is quite rapid, and it becomes faster for smaller values of 
the coupling constant. In particular, in the ^5*0 channel of the scalar meson exchange kernel 
it is sufficient to take into account Nmax ~ 4 — 5 terms in ©-(111) for the values of meson 
mass /i > 0.1 GeV and coupling constant < 40. In the direct calculations it also appears, 
that the Gaussian mesh with No = 64 is almost always enough, since an increase of No 
beyond this value practically does not change the results. 

Nevertheless, it should be pointed out that at certain conditions convergence of the 
solutions becomes poor or even breaks. In particular, it is lost at small meson masses ~ 0, 
and similar behaviour was found in [11] for the bound states of scalar particles. In this case, 
the introduction of form factors allows one to improve the situation. However, in total, we 
found that in general the equation ([1]) has stable solutions only for coupling constants 
below some critical value g"^^, which depends on the type of interaction and the channel 
considered. It can be found from the numerical calculations, since at coupling constants 
above some critical value the solution disappears, i.e. it becomes strongly dependent on 
Nmax, Ng and other numerical parameters. 

To find the critical value of the coupling constant the dependence of the ground state 
mass on the cut-off parameter A at fixed g^ have been investigated. The obtained results 
are presented in Figs. [3] and S] for ^Sq and ^Si-^Di channels respectively. It is evident 
that in the limit A — oo and coupling constants g^ < 40 the ^5*0 bound state mass does 
not depend on A. The same is valid for solutions (vertex functions). More exactly, the 
critical constant g'^^. for ^5*0 state is found to be 40.3. Theoretical estimation of the critical 
coupling constant can also be performed, and it gives the value g'^^. = 47r^ ^ 40. As it is 
seen from Fig. HJ the similar situation holds for the ^Si-^Di channel, where the critical 
value is g'^^. ~ 65. 

It can also be shown that g^ is directly connected to the power of decrease of vertex 
functions at large p, which is shown at Fig. [2l The solid line there corresponds to the 
solution below gcr with finite normalization, and the dashed one reproduces the solution 
beyond critical coupling, which seems to have infinite normalization. 
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4 Conclusion and acknowledgements 



A new method of solving the BS equations for the bound states of spinor particles by 
using the expansion of the vertex functions over the complete set of four-dimensional hy- 
perspherical harmonics is suggested. Within this method the BS equation is treated in a 
ladder approximation for the cases of scalar, pseudoscalar and vector meson exchanges with 
corresponding form factors. This method is shown to be effective and stable. 
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